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Abstract 

Macroscopic  viscoplastic  constitutive  models  for  y-y’  Ni-base  superalloys  typically  do  not 
contain  an  explicit  dependence  on  the  underlying  microstructure.  Microstructure  dependent 
models  are  of  interest  since  the  sizes,  volume  fractions,  and  morphologies  of  primary,  secondary, 
and  tertiary  precipitates  can  substantially  affect  the  stress-strain  response.  I  he  principle 
microstructural  features  that  can  significant!  ^ffect  the  stress-strain  response  of  y-y'  Ni-base 
supcralloys  are  the  grain  size  and  precipitate  volume  fraction  and  size  distributions.  An 
Artificial  Neural  Network  (ANN)  is  used  to  correlate  the  material  parameters  in  an  internal  state 
variable  cyclic  viscoplaslicity  model  with  these  microstructure  features  using  a  combination  of 
limited  experiments  augmented  by  polycrystal  plasticity  calculations  performed  on  other 
microstructures  within  the  range  characterized  experimentally.  The  trained  model  is  applied  to  an 
example  of  component  notch  root  analyses  to  explore  the  potential  impact  of  microstructurc- 
sensitive  constitutive  models  in  fatigue  design  of  structures. 
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1.  Introduction 


Conventional  macroscopic  (referred  to  as  ‘macro’)  models  for  cyclic  viscoplastic 
deformation  do  not  typically  contain  explicit  dependence  on  the  underlying  microstructure. 
Microstructure-sensitive  models  are  of  interest  in  view  of  variations  of  microstructure  within 
heat-treated  components.  Moreover,  components  with  functionally  graded  microstructures  may 
exploit  microstructure-property  relations  to  achieve  enhanced  fatigue  resistance.  A 
microstructure-sensitive  macroscopic  internal  state  variable  model  (McDowell  2005a,  2005b)  is 
useful  since  it  is  capable  of  giving  microstructure-sensitive  estimates  of  notch  root  cyclic  stress- 
strain  behavior  of  components. 

The  hierarchical  methodology  shown  in  Figure  1  represents  how  this  microstructural 
dependence  has  been  embedded  in  the  ’macro'  model.  In  this  work,  constitutive  models  are 
formulated  at  two  length  scales:  a  lower  scale  internal  stale  variable  (ISV)  dependent  microscale 
model,  which  is  formulated  at  the  length  scale  of  grains  using  the  crystal  plasticity  framework 
(cf.  Wang  et  al.  2007:  Shenoy  et  al.  2007).  and  a  macroscale  ISV  model.  Virtual  microstructure 
simulations  from  the  microscale  model  are  used  to  inform  the  microstructure  dependent 
parameters  in  the  macroscale  model.  There  are  three  primary  steps  in  the  algorithm  to  embed 
microstructure  dependence  into  the  ‘macro’  model  (Shenoy  2006;  Shenoy  et  al.  2007): 

•  Step  1  -  A  stress-strain  database  was  generated  for  different  strain  histories  using  the 
polycry  stal  plasticity  (referred  to  as  'micro')  model  for  a  range  of  microstructures  that 
are  intermediate  to  7-8  selected  actual  microstructures  used  for  ‘micro’  model 
calibration. 

•  Step  2  -  The  Walker-type  macroscopic  (‘macro')  model  parameters  were  determined 
for  the  stress-strain  data  for  each  microstructure  in  the  database,  including  virtual 
microstructures  intermediate  to  those  experimentally  characterized,  using  an 
optimization  scheme. 
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•  Step  3  -  An  Artificial  Neural  Network  (ANN)  was  trained  to  relate  the  parameters  in 
the  •macro’  model  to  the  corresponding  microstructural  parameters.  The 
microstructural  features  were  varied  randomly  within  the  range  shown  in  Table  1.  A 
total  of  150  microstructures  were  generated  to  train  the  ANN;  more  simulations  could 
be  run  if  desired. 

The  primary  microstructural  features  that  can  significantly  affect  the  cyclic  stress-strain 
response  of  Ni-hase  superalloys  are  the  grain  size  and  the  sizes  and  volume  fractions  of  the  Ni3AI 
y'  precipitates,  with  ranges  in  the  present  Ni-basc  superalloy  listed  in  fable  1.  Both  subsolvus 
and  supersolvus  heat  treatments  are  considered.  The  response  ot  this  particular  alloys  is  similar 
to  rN  100  (Milligan  el  al.  2004),  which  was  used  in  our  prior  study  (Shenoy  2006;  Shenoy  et  al. 
2007). 


Macroscale 
Parameter  Set 
(Ai  etc.) 


Microstructure 
Sensitive 
Macroscale 
Walker  type 


Parameters 


Parametric  Microstructure  Range 


Calibrated  'micro 
3-D  FE  model 


Material  Point 
Simulator  in  FE 
Code 


Parameter  Set  (f,,  d2,  f2,  d3,  dgraln) 


EPOGV  / 
Fitting 


Distributions  of  slip,  cumulative  plastic  strain 
residual  stress  among  grains  _ 


Material  Stress-Strain 
Response  input  into 
Lifing  Methodology 


Advanced  fatigue  life  models 
for  mean  and  variance 


Calibrated  'macro' 
Walker  -  type  ISV  model 


Figure  1.  Framework  for  linking  the  ‘micro’  model  to  the  ’macro’  model. 
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Table  1.  Ranges  of  microstructural  features  used  to  generate  the  database  for  150 

microstructures. 


Microstructural  Attribute 

Min 

Max 

Volume  Fraction  of  Primary  y'  Precipitate,  fpj 

0/0.07* 

0.25 

Volume  Fraction  of  Secondary'  ’/'  Precipitate,  fpo 

0.29 

0.46 

Secondary  y'  Precipitate  size,  d?  (nm) 

1 10 

340 

Volume  Fraction  of  Tertiary'  y'  Precipitate,  fp3 

0.024 

0.14 

Tertiary  y'  Precipitate  size,  d3  (nm) 

7.0 

21.0 

Grain  Size.  dgrajn  (microns) 

4.0 

34.0 

*0  pertains  to  the  case  of  no  y'  primary  precipitates:  otherwise,  minimum  is  0.07. 


2.  ‘Micro’  Polvcrystal  Plasticity  Model 

Local  crystal  plasticity  theory  is  used  at  the  grain  scale,  with  slip  accounted  for  on  both 
octahedral  and  cube  slip  planes  in  this  class  of  alloys.  The  kinematical  foundation  for  single 
crystal  elastoplasticitv  has  been  built  on  a  continuum  perspective.  (Asaro  1983a;  1983b).  The 
deformation  of  a  crystal  is  composed  of  contributions  from  (i)  an  overall  •“elastic”  distortion  of 
the  lattice,  and  (ii)  plastic  deformation  associated  with  dislocation  glide  that  does  not  disturb  the 
lattice  geometry.  The  deformation  gradient.  F.  is  multiplicatively  decomposed.  i.e.f 

F  -  F'»FP .  The  linear  hyperelastic  relation  (small  elastic  strains  are  assumed)  is  assumed  at 
each  temperature.  The  kinetics  of  shearing  on  each  slip  system  is  assigned  as  a  function  of 
resolved  shear  stress,  with  effects  of  initial  yield  and  work  hardening  incorporated  through 
micromechanical  relations  that  reflect  si/e-dependent  precipitate  shearing  versus  looping,  and 
dislocation  density  evolution.  The  model  equations  are  summarized  in  Table  2.  The  model  is 
integrated  over  a  large  number  of  grains  with  random  crystallite  orientation  distribution  function 
that  constitutes  a  Representative  Volume  Element  for  purposes  of  fitting  complex  cyclic 
deformation  behavior  of  actual  microstructures.  Details  may  be  found  elsewhere  and  are  beyond 
the  scope  of  this  paper  (Shenoy  2006:  Shenoy  et  al.  2007). 


-4- 


Table  2.  Equations  of  the  microstructure-sensitive  ('micro')  crystal  plasticity  model. 


Flow  Rule  with  Back  Stress  x“ .  Threshold  Stress  k“  and  Drau  Stress  Du 


r-  +r,  pi)  *»(*■-*■) 


where  the  threshold  stress  depends  on  dislocation  density  and  particle 
shear  resistance,  i.e., 

<  =  <>.  + «,  for  =  o^  cub 

<-  =  f  ( ■ < .  f  +  ( f„  •  d, .  rp! .  d,.f„  )■•  T'"‘  +  <, 

*“  ^  *p_rif 

r  n  n 

Ko.cub  =  I  (  To.  T  )  +  M,cub(fpl’^2»^p2’^3^p3)  J 
<u  =  hpcT^  +  hcb|x"|  +  h,cx“ 

V,*  =  V«*  =  cp, .  0.8 1  w  f;,  ny.  J-fL  +  cp2  Jo. 8 1  w  fp2  nr.  + 


f  ,\PB 


c  ,  III  -2^-f'2_aa.+c  d 

py  rd}nt  p3  b  *  * 

f  f  f 

f '  _  pi  fl  _  P2  f'  _  pj 

pl  fp,+fm’p2  fp2+f,/p3  f„+f. 


Internal  State  Variables  at  Slip  System  Level 
a)  Dislocation  Density 

P  r  =  I h<dW  {z0  +  k,  k2p^}|Y(f5)| 


zo  =  ^- 

hdjejr 


1  1 

OficIT  ~  .  +  , 

\Oia  Q:s 


h“p  =  l 


b)  Back  Stress 


xi"’=Cl(ilH.b>/pFsgn(tl-'-xrl)-xi“’]r<°1  .  1  = - ^ 

t  >  7.  +k.J 


Zo+Mp^ 
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3.  ‘Macro’  Cyclic  Viscoplasticity  Model 

3. 1 .  'Macro  ”  Model  Constitutive  Equations 

In  this  section,  the  focus  is  on  the  general  3D  formulation  of  the  unified  creep-plasticity 
internal  state  variable  macroscopic  model  used  to  calculate  the  stress-strain  response  for  different 
loading  conditions.  Under  isothermal  conditions,  the  rate  of  deformation  tensor  (true  strain  rate) 
D  is  decomposed  into  elastic  and  inelastic  components,  respectively  identified  by  the 

superscripts  e  and  n,  i.e..  D  =  D'  +  D” .  Under  small  strain  conditions,  the  stress  response  is 
given  by  the  material  time  derivative  of  the  Cauchy  (true)  stress,  i.e.,  q  =  2/rD''  +  ADrnI .  where 
l.i  and  a  are  the  Lame  constants.  Employing  a  two-term  potential  flow  rule,  the  inelastic  strain 
rate  is  w  ritten  as 


Dn  =  / 


with  a  two  term  flow  potential,  i.e.. 


CD 


(2) 


and 

5 = i?'  ■  «i = [(?' -«) : =  (?'- %)T  (3) 

Here,  q'  is  the  deviatoric  part  of  the  stress  tensor,  a  is  the  back  stress  tensor  (dcviatoric). 
and  q'-a  is  the  overstress.  The  exponents  N  and  Nz  indicate  the  rate  sensitivity  of  each  term 
in  the  flow  rule.  The  use  of  two  terms  in  the  flow  rule  facilitates  a  multi-mechanism  treatment  of 
plastic  flow  at  low  and  high  stress  levels.  We  set  N  >  N2  so  that  the  ilrst  term  in  the  flow  rule  is 
of  quasi  rate-independent  character  to  capture  the  dominant  cyclic  behavior;  the  threshold  stress 
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R  plays  a  role  in  this  dominantly  athermal  term.  The  second  term  largely  represents  effects  of 
thermally  activated  flow  at  low  stress  levels  and  higher  temperatures,  and  is  used  to  model  stress 
relaxation  behavior.  The  two  term  potential  may  in  itself  be  regarded  as  an  approximation  of  the 
power  law  breakdown  regime  at  higher  strain  rates  (cf.  Chaboche.  1989),  thereby  rendering  the 
exponential  terms  nonessential  (hence,  we  set  Bx  =  B{  =0 ).  The  constant  D  represents  drag 

stress,  and  k  is  the  threshold  stress.  The  material  hardening  derives  from  the  dynamic 
equilibrium  between  competing  hardening  and  dynamic  recovery  processes. 

The  term  R  corresponds  to  the  isotropic  hardening  of  the  material  (cf.  Chaboche.  1985).  i.e. 

R  =  Qh  i-cxp  ~bt,JjP  j-  P  =  JIMk'  (4) 

which  can  be  written  in  differential  form  as 

/?  =  (£,-/</?)  E,  =Qhbk,Fx  -bh  (5) 

where  Oh  is  the  saturation  limit  of  the  isotropic  hardening.  bh  controls  the  evolution  rate,  and  p 
is  the  normed  cumulative  inelastic  strain. 

Assuming  that  the  back  stress  tensor,  a ,  is  decomposed  into  two  components  (Moosbrugger 
and  McDowell  1989.  1990),  the  kinematic  hardening  relation  can  be  written  as 


where  A,  and  Bx  are  respectively  the  direct  hardening  and  dynamic  recover  rate  coefficients,  and 
a\m  and  is  the  saturation  limit  of  the  i,k  component  of  back  stress.  It  is  noted  that  more  complex 
forms  involving  nonlinear  dynamic  recovery  can  be  introduced  (cf.  Ohno  and  Wang  1991a. 


-7- 


1991b;  McDowell  1994,  1995;  Abdel-Karim  and  Ohno  2000)  if  necessary  to  model  cyclic 
ratcheting  or  progressive  mean  stress  relaxation  effects  more  accurately. 


3.2.  Finite  Element  Implementation  of  the  ‘Macro  ”  Model 

Either  explicit  or  implicit  integration  algorithms  can  be  used  for  time  integration  of  the 
foregoing  relations.  For  robustness  with  regard  to  time  step,  an  implicit  algorithm  procedure  is 
used.  A  residual  for  the  flow  rule  in  (l)-(3)  may  be  rewritten  as 


g,  =  4P,  ~ 


(7) 


The  residual  in  Equation  (7)  should  vanish  at  the  end  of  the  time  step.  To  do  this,  the  internal 


state  variables  (a,  R)  have  to  be  obtained  at  the  end  of  the  time  step,  which  in  turn  depend  on 


the  inelastic  strain  increment;  this  necessitates  an  iterative  procedure  is  needed  here.  A  Newton- 
Raphson  (N-R)  scheme  is  employed  for  the  root-finding  procedure,  w  ith  the  algorithm  described 
in  the  Appendix. 


4.  Microstructurcs 

A  matrix  of  potential  /  -/  microstructures  that  exhibit  different  average  grain  sizes  as  well 
as  different  sizes  and  volume  fractions  of  the  y‘  precipitates.  These  microstructures  are  listed  in 
Table  3.  The  definitions  of  the  microstructural  terms  are  listed  in  the  first  row  in  Table  3;  fpi,  fp:, 
and  fP3,  are  the  volume  fractions  of  the  primary',  secondary  and  tertiary  /  precipitates, 
respectively;  di.  d:,  and  dj  are  the  average  diameters  of  the  primary,  secondary  and  tertiary  / 
precipitates,  respectively,  and  dgMm  is  the  average  grain  size.  The  elastic  behavior  is  assumed 
linear  and  isotropic,  with  Young’s  modulus,  £mod  and  elastic  Poisson's  ratio,  v . 
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For  each  microstructure,  the  ‘macro’  model  parameters  were  obtained  using  the  trained  ANN 
(Shenoy  2006;  Shenoy  ct  al.  2007)  and  are  tabulated  in  Table  4.  The  Neural  Network  Toolbox 
in  MATLAB  was  used  to  perform  the  training  using  150  simulation  datasets.  In  this  particular 
ANN,  only  two  layers  were  used,  a  hidden  layer  and  an  output  layer.  The  microstructure 
parameters  are  used  as  the  input  parameters  and  the  ‘macro’  model  parameters  are  used  as  the 
output  variables. 


Table  3.  Features  of  model  microstructures. 


Microstructure  fpi 

di  (pm)  fp2 

d2  (nm) 

fp3 

d^  (nm) 

^grain 

(pm) 

M.1 

0.1 

1.0  0.45 

150 

0.1 

15 

31.8 

M.2 

0.1 

1.0  0.45 

250 

0.1 

15 

31.8 

M.3 

0.1 

1.0  0.3 

200 

0.2 

15 

31.8 

M.4 

0.1 

1.0  0.45 

250 

0.1 

15 

4.0 

M.5 

0.1 

1.0  0.45 

250 

0.1 

15 

16.0 

M.6 

0.1 

1.0  0.45 

250 

0.1 

15 

30.0 

Table  4. 

Macroscopic  viscoplastic  model  parameters  obtained  from  ANN  (650°C). 

v  4 

4 

Micro.  N 

(MPa)  (MPa) 

(MPa) 

(MPa) 

M.1 

95 

676  228000 

1041 

491 

6.5 

0 

95 

633  240000 

1114 

467 

6.5 

0 

M.3 

95 

772  208000 

960.6 

575 

6.5 

0 

M.4 

95 

755  278000 

1064 

499 

6.5 

0 

M.5 

95 

721  257000 

1114 

493 

6.5 

0 

M.6 

95 

646  242000 

1114 

468 

6.5 

0 

Micro. 

4 

D  4iwii 

V 

4 

4 

4  xIO'15 

(s') 

(MPa)  (MPa) 

(MPa) 

(s') 

M.1 

1.25  x  10‘J 

250  158800  0.33 

2256 

5.6 

132900 

7 

M.2 

1.25  x  1 0’3 

250  158800  0.33 

1640 

4.1 

210400 

7 
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M.3 

1.25  x  1 0'3 

250 

156400 

0.33 

2584 

6.5 

287700 

7 

M.4 

1.25x1  O'3 

250 

158700 

0.33 

1915 

4.8 

320500 

7 

M.5 

1.25  x  10’3 

250 

158700 

0.33 

1739 

4.3 

273700 

7 

M.6 

1.25  x  1 0'3 

250 

158800 

0.33 

1643 

4.4 

216800 

7 

5.  Geometry  and  Simulations 

5. 1  Notched  Component  Geometry.  Mesh  and  Element  Selection 

The  notched  component  (specimen)  used  for  all  the  simulations  is  specitled  in  Figure  2. 


Conlar  Line 


60* 


Figure  2.  Specified  geometry  of  notched  specimens  for  'macro'  model  FE  simulations,  show  ing 
details  of  notch  root  geometry  in  expanded  view  at  bottom. 
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Axisymmetric  elements  were  selected  for  the  analyses  as  they  provide  the  modeling  of 
bodies  of  revolution  under  axially  symmetric  loading  conditions.  Both  4-node  and  3-node 
continuum  axisymmetric  ABAQUS  elements  CAX4  and  CAX3,  respectively,  were  utilized  for 
all  of  the  simulations. 


Table  5.  Comparison  of  different  meshes  tested  to  ensure  convergence  of  the  solution  at  the 

notch  root. 


Mesh 

Elements 

Axial 

Component 

of  Plastic 

Strain  (After 

10  cycles) 

Element 

Volume  at 

Notch  Root 

(pmJ) 

Stress 

Concentration 

Factor 

(Elastic) 

a 

314 

1.00% 

9.9 

b 

421 

1.15% 

38000 

10.5 

c 

556 

1.27% 

16900 

10.9 

d 

665 

1.30% 

12400 

11.1 

e 

789 

1.34% 

9500 

11.2 

f 

1045 

1.37% 

6000 

11.5 

9 

1399 

1.40% 

4200 

11.6 

h 

2927 

1.45% 

1500 

11.9 

i 

5750 

1.47% 

380 

12.1 

1 

3178 

1.47% 

200 

12.1 

Final 

1055 

1.46% 

670 

11.7 

The  component  defined  in  Figure  2  was  modeled  as  axisymmetric.  Several  different  meshes 
were  compared  to  ensure  that  accurate  results  were  obtained  at  the  notch  root.  As  seen  in  Table 
5.  the  results  converged  when  the  area  of  the  elements  at  the  notch  root  were  approximately  400 
pm2.  It  is  noted  that  even  though  the  area  of  the  elements  at  the  notch  is  less  for  each  successive 
mesh  from  mesh  a  to  j  (see  Tabic  5).  the  number  of  elements  were  reduced  between  meshes  i 


and  j  by  coarsening  elements  further  away  from  the  notch  to  reduce  computation  lime.  In  order 
to  reduce  computation  time  further,  a  mesh  was  selected  with  only  1055  elements  using  similar 
methods  of  fanning  out  the  mesh  across  less  critical  regions:  however,  the  area  of  the  elements  at 
the  notch  root  were  kept  such  that  the  maximum  axial  component  of  the  plastic  sttain  was  within 
1%  of  the  result  of  the  finest  mesh  tested.  This  mesh  can  be  seen  in  Figure  3.  The  stress 
concentration  factors  as  determined  by  the  FEA  simulations  in  the  elastic  regime  are  also  shown 
in  Table  5.  An  analytical  elastic  solution  for  the  supplied  geometry  was  not  available. 


Figure  3.  The  final  selected  mesh  used  to  model  the  notched  component:  a)  complete  mesh,  and 

b)  mesh  in  the  vicinity  of  the  notch  root 


5.2  Finite  Element  Simulations 

FE  Simulations  w'ere  performed  using  the  ’macro’  UMAT  for  all  the  microstructures.  For  all 
simulations,  isothermal  conditions  were  assumed  at  650°C.  In  addition,  the  FE  simulations  were 
all  performed  under  nominal  (in  grips)  stress  control  with  R  -  0.05.  The  load  histories  used  in 
the  FE  simulations  are  listed  in  Table  5.  Notch  root  stress-strain  responses  and  associated 
accumulated  plastic  strain  for  all  the  microstructures  arc  given  in  Figures  4-1 1.  The  accumulated 
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plastic  strain  at  each  point  is  defined  as  j dt ,  with  e!  =  tyDy  .  Calculations  of  nonlocal 

measures  of  notch  root  stress  and  strain  were  performed  in  all  cases  by  averaging  the  values  over 
a  small  but  finite  area  of  approximately  2500  pm2  to  provide  an  indication  of  the  driving  force 
for  formation  of  small  fatigue  cracks  of  finite  size,  and  to  mitigate  mesh  dependencies  of  the 
solustions.  These  are  the  stress  and  strain  values  reported  in  Figures  4-1 1. 

The  sensitivity  of  responses  to  mierostructure  variation  is  interesting.  For  the  shorter  dwell 
times  the  creep  resistance  of  microstructures  (i.e..  lowest  accumulated  plastic  strain),  listed  in 
decreasing  order,  were  M.4.  M.3,  M.5,  VI.  1.  V1.6.  and  M.2.  respectively.  However  for  the  longer 
hold  times,  the  microstructures  with  the  best  creep  resistance  were  VI.  1,  M.4.  M.6.  VI. 2.  M.5.  and 
VI. 3,  respectively.  The  largest  difference  in  behavior  between  short  and  long  hold  times  was 
observed  for  microstructure  M.l,  where  the  smaller  size  of  the  secondary  precipitates  and  the 
larger  grain  size  yielded  better  creep  behavior  for  longer  hold  times  than  at  shorter  hold  times. 
Microstructure  M.l  had  lower  accumulated  plastic  strain  for  a  short  hold  time. 


Table  6.  Strain  histories  used  in  FE  simulations  to  compare  response  of  various  microstructures. 


History 

Temp 

(F) 

Applied 

Nominal 

Max. 

Stress 

(MPa) 

R 

Ramp 

Up  (s) 

Ramp 

Down 

(s) 

Peak 

Dwell 

(s) 

Base 

Dwell 

(s) 

CPM 

1 

1200.0 

200.0 

0.05 

1.5 

1.5 

0.0 

0.0 

20.000 

2 

1200.0 

200.0 

0.05 

15.0 

15.0 

90.0 

0.0 

0.500 

3 

1200.0 

200.0 

0.05 

15.0 

15.0 

1200.0 

0.0 

0.049 

4 

1200.0 

200.0 

0.05 

15.0 

15.0 

12000.0 

0.0 

0.005 

13  - 


_  2.5 
£ 

S 

<o 
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Time  (s) 

Figure  5.  Comparison  of  notch  root  accumulated  plastic  strain  (  e '  dt )  of  all  tested 
microstructures  for  microstructures  M.I-M.6  for  History  1 . 
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Figure  6.  Comparison  of  notch  root  axial  true  stress  ( er13 )  versus  axial  true  strain  ( eu  )  of  all 
tested  microstruelures  for  microstructures  M.1-M.6  for  History  2. 


Figure  7.  Comparison  of  notch  root  accumulated  plastic  strain  (  [  t '  df )  of  all  tested 
microstructures  for  microstructures  M.1-M.6  for  History  2. 
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8000  10000  12000  14000 


Figure  9.  Comparison  of  notch  root  accumulated  plastic  strain  (  f  eHdt )  of  all  tested 


microstructures  for  microstructures  M.1-M.6  for  History  3. 
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Figure  10.  Comparison  of  notch  root  axial  true  stress  ( crn  )  versus  axial  true  strain  (e22 )  of  all 
tested  microstructures  for  microstructures  M.1-M.6  for  History  4. 


Figure  1 1.  Comparison  of  notch  root  accumulated  plastic  strain  (  £  c  dt )  of  all  tested 
microstructures  for  microstructures  M.1-M.6  for  Hstory  4. 
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Conclusions 


Two  important  observations  were  made  from  this  study: 

•  Using  an  artificial  neural  network  to  determine  'macro'  material  parameters  to 
facilitate  simulations  of  various  microstructures  is  potentially  useful  in  exploring 
effects  of  possible  microstructure  compositions  on  service  behavior  of  components. 
The  response  of  such  ‘engineered’  microstructures  can  be  projected  before  such 
materials  are  physically  realized.  Such  studies  can  enable  setting  target  volume 
fractions  and  precipitate  sizes  in  the  case  of  Ni-base  superalloys. 

•  The  sensitivity  study  of  variation  of  notch  root  response  with  microstructures  M.l- 
M.6  provided  insight  into  the  potential  of  the  method  outlined  here  for  linking 
polycrystal  plasticity  calculations  to  macroscopic  constitutive  relations  for  use  in 
notch  root  analyses.  There  is  a  pronounced  effect  of  microstructure  on  mean  stress 
relaxation,  for  example,  that  may  impact  fatigue  life. 
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APPENDIX  -  Details  of  Numerical  Implementation 


Starting  with  the  residual  in  Equation  (7),  we  frame  the  N-R  iteration  as  follows: 

g, 


4P/+,  «  A/;, 


[dg/dAp\ 


The  flow  rule  can  be  rewritten  as 


=  4P"4  (q'.q.R)- AVj  ^AtF:  (q'.q) 


where  F(q‘, a,R)  =  ^  • 

'  '  V  / 


From  chain  rule  differentiation. 


dg_ 

5A/? 


,  ,  !J  f  8F  da'  dF  da  dF  dR  "l 

=  l-A.-Af  : — — + — : — — + - 

'  V  2  dq  dAp  da  dAp  SR  dAp  J 

•'  V  2  dq  dAp  dq  dAp  J 

(  2  fj  ^2 

2pn  +  j(Al  +  A2)n-]jj(Rlql  +  B2q2)  +j(E,-F,)R 

2  [2 

2  /in  +  j(4  +  AI),j-jj(Blql  +  B2q2 ) 


=  1  +  F 


n : 

L~  V 


where 


’  2  D  [  D 


N-l  f 

„  3  .  .  A, 

F  =  —  A  At  — 

‘•2  2  ;  D 


m 

D 


,V,-I 


(Al) 


(A2) 


(A3) 


(A4) 


The  inelastic  strain  increment  and  the  internal  state  variables  arc  updated  using  this  scheme 
through  successive  iterates  until  gM  ~  0  at  the  end  of  the  time  step. 

The  Newton-Raphson  technique  achieves  a  rapid  convergence  rate  near  the  actual  solution. 
However,  it  may  diverge  if  the  initial  guess  is  sufficiently  far  from  the  actual  solution.  This 


-  19- 


problem  occurs  when  the  values  of  D n  change  significantly  over  a  time  step  due  to  a  large  time 
step.  Accordingly,  a  line  search  method  proposed  by  Cuitino  and  Ortiz  (1992)  is  employed. 

In  addition,  a  time  step  subincrementation  algorithm  can  also  be  supplemented  to  ensure 
convergence  based  on  the  work  by  Bennett  (1999)  and  McCiinty  (2001).  The  implementation  of 
the  subincrementation  method  consists  of  •‘wrapping"  a  lime  step  subincrementation  routine 
around  the  Newton-Raphson  and  line  search  methods.  If  the  Newton-Raphson  with  line  search  is 
not  satisfactorily  converging  for  a  given  time  step,  then  the  integration  is  performed  in  two 
smaller  time  steps,  obtained  by  halving  the  original  time  step.  The  flow  chart  of  this  procedure  is 
given  in  Figure  12. 

The  viscoplastic  'macro'  constitutive  model  is  coded  as  an  ABAQUS  6.5-4  User  MATerial 
subroutine  (UMAT).  ABAQUS  provides  the  strain  and  Jacobian  (the  derivative  of  the  stress  with 
respect  to  elastic  strain)  matrices  at  the  start  of  the  time  step,  along  with  the  overall  strain  and 
time  increments  for  the  step.  The  stresses.  Jacobian  matrix,  and  internal  variables,  which  are 
history  dependent,  must  be  updated  at  the  end  of  the  time  step  in  this  UMAT.  These  values  are 
returned  to  ABAQUS,  which  then  continues  equilibrium  iteration.  It  is  noted  that  the  back 
stresses  are  not  rotated  at  the  end  of  the  time  step  as  would  be  required  in  a  finite  strain 
implementation,  owing  to  relatively  small  cyclic  strains  of  perceived  applications.  This  could  be 
added  if  necessary. 
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Begin  with  overall  strain  and 
time  steps.  Afr?  A/ 


Figure  12.  Flowchart  of  the  constitutive  law  update  in  ABAQUS  UMAT  . 


It  is  necessary  to  formulate  the  Jacobian  matrix  required  by  ABAQUS  at  the  end  of  the  time 
step  for  global  equilibrium  assessment.  The  Jacobian  J  is  given  by 


-21  - 


J  = 


dAa  8 


dAe  dAe 


[C:(£>-D”)A/]  =  C-(C:»)<8> 


dAp 

dAe 


(A5) 


The  magnitude  of  the  increment  of  plastic  flow  can  be  written  as 


Ap  =  A  yf^AtF(<T',q,  R)+A,t  ^|a/F:  {q'.q) 


(A  6) 


Chain  rule  differentiation  gives 
dAp 


dAe 


.  [J.fdF  da'  8F  da  ~dAp  dF  dR  ~dAp 
\  2  ^  da  dAe  8q  dAp  dAe  dR  5Ap  dAe 


dF:  dq  |  8FZ  dq  dAp' 
dq  dAe  dq  dAp  dAe  / 


(A7) 


The  contribution  of  the  deviatoric  stress-strain  response  to  the  Jacobian  can  be  assessed  by 


a'  =  a  —  aiiI  => 


8Aq‘  _  dAq  d 
dAe  dAe  dAe 


(K&sul_) 


(AS) 


Since  the  last  term  is  purely  elastic  response. 


^-  =  C-(C:n)(S>^—  K(I®I) 
dAe  ~  v-  dAe  u 


where  C  is  the  elasticity  tensor.  Gathering  results  leads  to  the  Jacobian 


(A9) 


7=C-(C:«)0 


I  +  F,  +  F, 


(A10) 


where 


F,  =  K 


Fa=J\. 


"  :|c:,  +  ^1  +  4) 2  “  ) 


(All) 
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